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Positronium hyperfine splitting: analytical value at 0(ma 6 ) 
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We present an analytic calculation of the 0(raa 6 ) recoil corrections to 
the hyperfine splitting (HFS) of the ground state energy levels in positro- 



nium. We find AE iec = ma 6 (-± In a + f§ - ^ - ^ + ^) 



ma 6 (— | In a + 0.37632 J , confirming Pachucki's numerical result We 



present a complete analytic formula for the 0(ma & ) HFS of the 

positronium ground state and, including 0(moJ In 2 a) effects, find 

(£(l 3 Si) - ^PSo)) thcory = 203 392(1) MHz. This differs from the experi- 
mental results by about 3 standard deviations. 
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Spectroscopy of positronium, an atom consisting of an electron and a positron, provides 
a sensitive test of the Quantum Electrodynamics (QED) applied to bound state problems. 
Electron and positron are so much lighter than the lightest hadrons that the effects of strong 
interactions are negligible compared with the accuracy of present and any conceivable future 
experiments. For this reason positronium represents a unique system which can, in principle, 
be described with very high precision by means of the QED only. One should also mention 
that the measurements of positronium spectrum are performed with very high accuracy [F2J. 

There are two main approaches used in the studies of bound states. The Bethe-Salpeter 
method is based on an exact two-body relativistic wave equation . The other approach is 
the so-called Non- Relativistic Quantum Electrodynamics (NRQED) || , which is an effective 
field theory based on the QED for small energies and momenta. Thus, by construction, the 
NRQED takes advantage of non-relativistic energy of the electron and positron in positro- 
nium. 

We note in passing that similar techniques are used nowadays for describing heavy quark- 
antiquark bound states. From this perspective, positronium may serve as a testing ground 
for methods which can in the future be applied to the QCD. 

The HFS of the positronium ground state (i.e. the difference between the energies of 
the ground state with total spin 1 and 0) belongs to one of the most accurately measured 
physical quantities. Two experimental values of the highest precision are: 



obtained respectively in |5]||. The bulk of this effect is of the order ma A , where m = 
0.51099907(15) MeV (corresponding to 1.2355898(4) x 10 14 MHz) is the electron mass and 
a = 1/137.03599959(51) is the fine structure constant. Higher order corrections must be 
included to fully exploit the experimental accuracy. In particular, since ma 6 = 18.658 MHz, 
a complete calculation at this order is required. With an exception of the leading logarithm, 
effects of order ma 7 have not yet been studied. Clearly, the experimental precision warrants 
further studies of such corrections. 

The history of theoretical calculations of various contributions to the HFS of positronium 
is quite long. They can be represented by a series in powers and logarithms of the fine 
structure constant, 



Au = £(l 3 Si) 



E(VS ) = 203 387.5(1.6) MHz, 



(1) 



and 



Au 



203 389.10(0.74) MHz, 



(2) 




(3) 



The leading order O {ma 



4 ) HFS was obtained in [[/J-^. 
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(4) 




The first correction was calculated in [1C], 




(5) 
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The second correction consists of the following contributions: 



ma 6 n 2 = AE g _ 2 + A£ annih + A£ rad rec + AE rec . (6) 

The logarithmic contributions at this order, 0(ma 6 In a), present in the annihilation Ai? ann ih 
and recoil corrections AE rec , were found first [pT||12l . AE g _ 2 arises from the anomalous 
magnetic moment of the electron at 0(a,a 2 ). The three-, two- and one-photon annihila- 



tion contributions giving AE^^ were found in []13]-|15||, respectively. The non-annihilation 
radiative recoil contributions AE Ia d iec were studied in [ TO , [T7| , while pure recoil corrections 
AE rcc were discussed in 

For most of these contributions, several independent calculations were performed and 
an agreement was achieved. Moreover, the results for all contributions to HFS are known 
in the analytic form, with the exception of the pure recoil corrections AE rec . By pure 
recoil corrections one understands those induced by diagrams where each virtual photon is 
created by electron and absorbed by positron, as shown in Fig. [I} For these effects, three 
independent calculations arrived at three different results PU^ISf . The discrepancy has not 
been clarified so far and the resulting uncertainty in the theoretical prediction for the HFS 
of the ground state is much larger than the experimental error. The importance of clarifying 
this theoretical point has been emphasized by several authors. In this Letter we present an 
analytic calculation of these corrections. Numerically our result coincides with Ref. ]IJ. 

We start with a short description of the framework of our calculation, leaving the details 
to a separate publication. First, we calculate the on-shell scattering amplitude for non- 
relativistic (» < 1) particles to the necessary order. Along with the leading amplitude of 
a single Coulomb exchange, it includes the relative order 0(v 2 ) Breit corrections and also 
higher order 0(v 4 , otv 3 ) terms. By construction, it is gauge-invariant. Taken with a minus 
sign, this amplitude provides the effective potential for non-relativistic particles. 

Further, we solve the Schrodinger equation incorporating corrections to the Coulomb 
potential using ordinary quantum mechanical perturbation theory. According to standard 
rules, we get the 0(ma 6 ) correction to the ground-state energy as a sum of the relativistic 
corrections to the tree level and one-loop scattering amplitude, and of the second order 
correction due to the Breit potential. Previously, this scheme was used for the calculation 
of the 0(ma 6 In a) corrections to the levels of S-states |19[ and of the 0(ma 6 ) corrections 
to the levels of P-states [pOf| . 

An implementation of this program leads to a divergent result. The reason for this 
divergence is well known - it is the application of the non-relativistic expansion in the 
relativistic momentum region, where it is not appropriate. The divergence is canceled if one 
includes additional short-distance or hard-scale contributions to the scattering amplitude, 
which cannot be obtained from the non-relativistic expansion. 

To deal with the divergences in both the non-relativistic region and in the short-distance 
corrections we employ dimensional regularization. In the context of bound state calculations 
in QED this regularization scheme was used in Ref. ||21|| , where the known results for (D(ma 5 ) 
corrections to positronium energy levels were successfully reproduced. The advantage of 
the dimensional regularization is that it makes the matching calculation of the low-scale 
effective theory and the complete QED extremely simple. To obtain the contribution of a 
given Feynman diagram to the Wilson coefficient of the 5(r*)-like effective operator, we only 



3 



need to calculate that diagram for zero incoming momenta of all particles. We stress that 
this is only correct if one uses dimensional regularization for both infrared and ultraviolet 
divergences. With any other regularization scheme an additional calculation is required. 

We find that in the sum of the short- and long-distance contributions the singularities 
l/e = 2/(4 — D) disappear and one arrives at a finite result. 

Since the dimensional regularization is used throughout the Letter, we mention how the 
spinor algebra was treated. To calculate the shift in the ground state energy due to some 
operator Oi one has to calculate the trace of the form Tr , where $ is an appropriate 

wave function. The spinor parts of the relevant wave functions are: 

for para- and orthopositronium states, respectively. In the latter case £ is the polarization 
vector. The traces are calculated in the D-dimensional space. Since we always encounter 
an even number of 75's, we treat them as anticommuting. We also average over directions 
of the vector £. In order to obtain corrections to the HFS we first calculate separately the 
traces for ortho- and parapositronium states and then take the difference of the two. 

The problem is naturally divided up into the calculation of the matrix elements of the 
effective operators (soft contributions) and the Wilson coefficients of the effective 5(f)-like 
operators (hard contributions) in the effective Hamiltonian: 

^E rcc A so ft-E< rec -|- /^hard-^'rec- (7) 

The calculation of Wilson coefficients is always done for the incoming and outgoing particles 
at rest. Both technically and conceptually, this is close to the calculation of the match- 
ing coefficient of the vector quark-antiquark current in QCD and its NRQCD counterpart, 
described e.g. in p^j23|. 

This technique is remarkably useful for the so-called radiative recoil corrections to the 
HFS, where one of the three exchanged photons is created and absorbed by the same particle, 
as shown in Fig. ^|. It is sufficient to calculate the corresponding integrals exactly at the 
threshold in dimensional regularization, since there are no non-relativistic contributions to 
the radiative recoil corrections and no matching is required. Performing this calculation we 
obtain: 



A£ rad rec = ma 6 ^-4 + - In 2 + . (8) 

raarec | n 2 „ .„ „„ 2 v I 



'cw + ! ln2 _!£ + « 

2tx 2 3 48 36tt 2 

This result is in complete agreement with the analytic result published previously |TJ 

Applying the same technique to obtain the hard-scale contribution to the recoil correc- 
tions AE Tec we obtain :[] 



HQ? ( I 

M,ard£rec = ^l^ )! ( ~~ + 41nm - + ^ - 6 In 2 ) . (9) 




x We neglect factors T 2 (l + e) and (47r) 2<E which do not contribute to the final, finite result. 
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In the above equation ^rf(O) stands for the value at the origin of the ground state solution 
of the (^-dimensional Schrodinger equation. 

The calculation of the soft scale contributions requires the treatment of the relativistic 
corrections to the tree level and one-loop scattering amplitudes, as well as the second iter- 
ation of the Breit potential. The main difficulty associated with this calculation is that it 
should be done in d = 3 — 2e dimensions, thus necessarily spoiling some simplifying features 
of the Coulomb problem in three dimensions. Still, the calculation is feasible. Since the 
non-relativistic Hamiltonian is only singular for r — > 0, it turns out possible to extract this 
divergence in the form 1^(0) | 2 /e, without solving the Schrodinger equation in (i-dimensions. 
Our final result for all non-relativistic contributions reads 

A soft £ rec = ^l^(0)| 2 (~ - A\n{ma) + ^) . (10) 

In the sum of the hard and non-relativistic contributions, Eqs. (j9yi0"l), the 1/e divergences 
disappear and we can take the limit e — > in the sum. We thus arrive at the final result for 
the recoil corrections to the HFS of the positronium ground state: 

. _ . _ 6 / L 331 ln2 17C(3) 5 \ 

AE Tec = A hard £ rec + A soft £ rec = ma 6 [-- In a + — - — + -—r . 11 

\ 6 432 4 8ti z 127T J 

Numerically this is AE mc = ma 6 (— |lna + 0.37632^ which is in excellent agreement with 
Ref. IH, where for the non-logarithmic part of the correction a number 0.3767(17) was ob- 
tained. In view of the fact that in Ref. p| a different regularization was used, this agreement 
gives us confidence in the correctness of the result^. 

The recoil correction was the last correction to positronium bound state HFS not known 
analytically. Having obtained its value (Eq. (0)), we are now in position to present the 
final analytic result for the HFS of the positronium ground state including 0(ma 6 ) terms: 

E^S,) - E^So) = ma' - £ (j + Un2) 

« 2 r 5 2, 1367 5197 2 /221 2 1\ , „ 53^ /oN ll 

+ 7r 2 lna + tt 2 + tt 2 + - n2 C(3) }. (12) 

7T 2 [ 24 648 3456 V144 2J 32 sv; Jj v ; 

Numerically this corresponds to Au = 203 392.899 MHz, if we use the central values for m 
and a, given after Eq.(|^). 

To arrive at the final prediction for the HFS splitting of the positronium ground state, one 
should try to quantify the theoretical error. The errors caused by the uncertainties in the fine 
structure constant and the electron mass are ~ 0.003 and 0.07 MHz, respectively. The main 
uncertainty comes from the unknown higher order effects. Though formally ma 7 ~ 0.1 MHz, 
the leading 0(ma 7 ln 2 a) terms contribute —0.92 MHz to the HFS ||25|| . Therefore, it remains 



2 When this work was completed, we were informed |24|] about an independent numerical calcula- 
tion of the recoil corrections. Though that work is still in progress, its preliminary results seem to 
coincide with the results of Ref. Q and of the present work with rather good accuracy. 
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very important to calculate the remaining, non-leading terms in 0(ma 7 ). In this context 
we note that the complete 0(ma 6 ) correction, including the ma 6 In a; term, gives a shift 
of 11.79 MHz, whereas the term ma 6 In a alone contributes 19.12 MHz. We see that an 
estimate of the complete correction based on the ma 6 \n(a) approximation would not be 
accurate. 

At the moment the best we can do is to take the leading log contribution 0(ma 7 In 2 a) as 
an estimate of the higher orders corrections to HFS and thus attribute ~ 1 MHz uncertainty 
to the theoretical prediction. Combining this theoretical uncertainty with Eq. ([12]) , we 
obtain the theoretical estimate for the HFS of the ground state of the positronium: 

Az/ theory = 203 392(1) MHz. (13) 

Compared to the experimental results Eq. (^,@) we observe a significant deviation of the 
order of 3 — 4 experimental errors. We look forward to future improved measurements of 
positronium HFS and their confrontation with QED. 

We are grateful to A. Burichenko for informing us about his results prior to publication. 
K.M. would like to thank the High Energy Theory Group at the Brookhaven National 
Laboratory for hospitality during the final stage of this project. This research was supported 
in part by DOE under grant number DE-AC02-98CH10886, by BMBF under grant number 
BMBF-057KA92P, by Graduiertenkolleg "Teilchenphysik" at the University of Karlsruhe 
and by the Russian Foundation for Basic Research under grant number 98-02-17913. 
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FIGURES 




FIG. 1. Feynman diagrams representing pure recoil corrections to positronium HFS. 




FIG. 2. Examples of radiative recoil corrections to positronium HFS. 
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